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SPHERICAL-WAVE  CHARACTERIZATION  OF  INTERIOR 
AND  EXTERIOR  ELECTROMAGNETIC  SOURCES 


David  A.  Hill 

Electromagnetic  Fields  Division 
National  Institute  of  Standards  and  Technology 
Boulder,  Colorado  80303 


This  report  presents  spherical  scanning  formulations  for  the  general  case  where 
electromagnetic  sources  are  present  both  inside  and  outside  the  measurement 
region.  The  fields  are  expanded  in  terms  of  outgoing  waves  (due  to  interior 
sources)  and  standing  waves  (due  to  exterior  sources).  Additional  information  is 
required  to  solve  for  the  increased  number  of  unknowns,  and  results  are  derived 
for  the  spherical  wave  coefficients  in  terms  of  the  tangential  electric  and  magnetic 
fields  on  a sphere.  The  special  case  of  an  electrically  small  emitter  is  examined  in 
detail,  and  several  formulations  are  presented  for  the  components  of  electric  and 
magnetic  dipoles.  The  dipole  formulations  are  intended  for  use  in  detection  of  a 
weak  emitter  (such  as  a timing  device)  in  a noisy  environment. 

Key  words:  electric  dipole;  electrically  small  source;  magnetic  dipole;  near-field 
measurement;  spherical  scanning;  spherical  wave;  test  volume. 


1.  INTRODUCTION 

In  spherical  near-field  antenna  measurements  [ 1 ] or  in  characterization  of  electromagnetic 
interference  (EMI)  sources  [2],  the  radiated  fields  are  determined  outside  a minimum  sphere 
enclosing  the  sources.  For  radar  cross  section  (RCS)  or  receiving  antenna  measurements  [3,4], 
the  sources  are  located  outside  the  test  volume.  If  the  test  volume  is  spherical,  the  interior  fields 
can  be  expanded  in  spherical  waves.  In  the  case  of  either  interior  or  exterior  sources,  the  problem 
is  to  determine  the  coefficients  of  a spherical  expansion  from  measurements  made  with  a probe 
antenna  on  the  surface  of  a sphere. 

In  the  more  general  case,  sources  are  located  both  inside  and  outside  a spherical-shell, 
source-free  region.  This  situation  can  arise  in  practice  because  of  the  presence  of  unwanted 
sources  or  scatterers  located  either  inside  the  inner  sphere  boundary  or  outside  the  outer  sphere 
boundary  of  the  shell.  In  this  case,  the  problem  is  to  determine  the  coefficients  of  a spherical- 
wave  expansion  which  is  valid  in  the  spherical-shell  region.  More  information  is  required  because 
the  expansion  includes  both  outgoing  waves  (due  to  the  interior  sources)  and  standing  waves  (due 
to  exterior  sources).  The  purpose  of  this  report  is  to  present  the  formulation  for  this  more  general 
case  where  both  interior  and  exterior  sources  are  present.  For  typical  RCS  or  receiving  antenna 
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measurements  (where  the  intentional  source  is  outside  the  shell),  the  test  volume  may  be 
electrically  large,  and  the  interior  source  could  be  due  to  scattering  from  a support  structure 
located  inside  shell  region.  For  applications  to  characterizing  typical  EMI  sources  or  detecting 
weak  radiators  (such  as  timing  devices  for  explosive  detonators),  the  source  is  normally 
electrically  small,  and  the  spherical  wave  expansion  for  the  radiated  field  can  be  reduced  to 
electric  and  magnetic  dipole  terms  [2],  When  the  radiator  is  weak,  external  noise  sources  can 
compete  with  the  desired  signal,  and  they  can  be  modeled  as  the  exterior  sources  outside  the 
spherical  shell  region. 


2.  SPHERICAL- WAVE  EXPANSIONS 

In  source-free  regions  with  spherical  boundaries,  the  electric  and  magnetic  fields  can  be 
expanded  in  spherical  waves  which  satisfy  Maxwell's  equations  term  by  term.  Various  notations 
are  available  for  spherical  waves  [5,6],  but  this  report  follows  the  notation  in  Hansen  [1]  because 
he  has  a number  of  useful  mathematical  results  (orthogonalities,  radiated  power,  dipole  fields, 
etc.)  already  worked  out.  For  simplicity,  the  special  cases  of  interior  or  exterior  sources  only  are 
considered  in  sections  2.1  and  2.2.  The  general  case  of  interior  and  exterior  fields  covered  in 
section  2.3  can  be  considered  to  be  a superposition  of  the  two  special  cases. 

2. 1 Interior  Sources  Only 

This  section  treats  the  special  case  shown  in  figure  1 where  the  sources  are  confined  to  the 
interior  region,  r < r,.  In  the  source-free  region,  r ^ /*„  the  electric  field  E and  magnetic  field  H 
can  be  written  as  expansions  [1]  of  outgoing  waves: 
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The  time  dependence,  exp(-/a)/),  is  suppressed,  the  wavenumber  k = o)v/^€,  the  intrinsic 
admittance  r\  = y^e/ji,  and  e and  p are  the  permittivity  and  permeability  of  the  medium  (usually 
free  space).  Standard  spherical  coordinates  (/',6,4>)  are  used  throughout. 

The  transverse  electric  (TE)  spherical  wave  F^is  given  by  [1] 
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where  Pj"  is  the  normalized  associated  Legendre  function  [1,7]  and  is  a spherical  Bessel 
function  . For  c = 3,  the  spherical  Bessel  function  is 

- h^\kr),  (4) 


where  is  the  spherical  Hankel  function  of  the  first  kind  [8],  The  transverse  magnetic  (TM) 
spherical  wave  is  given  by  [1] 
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where 


The  unknown  coefficients  remain  to  be  determined.  In  theory  the  n summation  runs  from  1 
to  oo,  but  in  practice  the  summation  can  be  truncated  at  N somewhat  greater  than  kr^. 

The  spherical  waves  are  dimensionless,  and  the  normalization  is  chosen  so  that  the 
power  P radiated  through  the  sphere  of  radius  r,  (or  any  larger  sphere)  enclosing  the  source  is 
given  by  [1] 
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P = 


^ smn 


where  is  shorthand  for  the  triple  sums  in  eqs  (1)  and  (2).  The  units  of  are 

snot 

Spherical  near-field  scanning  methods  for  determining  the  coefficients  are  well 
known,  including  probe  correction  [1,3].  Probe  correction  is  not  considered  in  this  report,  and 
the  results  for  will  be  given  in  terms  of  field  components  (which  could  be  measured  with 
dipole  probes). 

In  order  to  determine  the  coefficients  from  the  tangential  components  of  the  electric 
or  magnetic  field,  an  orthogonality  relationship  between  the  tangential  components  of  the 
spherical  harmonics  is  needed.  Hansen  gives  such  a relationship,  eq  (A1.69)  in  [1],  but  it  is 
convenient  to  write  it  in  a more  compact  form.  The  first  step  is  to  define  tangential  spherical 
harmonics  that  are  independent  of  r: 
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Equations  (8)  and  (9)  are  related  by  the  following  cross  product: 
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Hansen's  orthogonality  relationship,  eq  (A1.69),  can  now  be  written  as 
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where  the  solid  angle  integration  over  4n  steradians  is  shorthand  for 
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Consider  now  the  case  where  the  tangential  electric  field,  = E^O  + E^$,  is  specified 
(measured)  over  the  surface  of  a sphere  of  radius  r (>  r,).  If  the  dot  product  of  a tangential 
spherical  harmonic  with  eq  (1)  is  integrated  over  angle,  eq  (1 1)  can  be  used  to  reduce  the 
summation  to  a single  term: 


/ / (13) 
An  yr) 

Since  R^(kr)is  either  a spherical  Hankel  function  (for  a = 1)  or  its  derivative  (for  a = 2),  it  has 
no  zeros  for  any  real  values  of  kr,  and  eq  (13)  can  be  solved  for  If  the  subscripts  (o,p,v) 

are  replaced  by  (s,m,n),  eq  (13)  can  be  rewritten  as 
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This  is  the  desired  result  for  the  coefficients  in  terms  of  the  tangential  electric  field,  and  it  is 
consistent  with  Hansen's  result  [1],  A similar  procedure  can  be  carried  out  starting  with  eq  (2)  for 
the  magnetic  field  to  obtain  the  coefficients  in  terms  of  the  tangential  magnetic  field; 
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The  tangential  electric  or  magnetic  fields  required  in  eqs  (14)  or  (15)  could  be  obtained  in 
theory  by  sampling  with  ideal  electric  or  magnetic  dipole  probes.  An  electrically  short  linear 
antenna  is  a good  approximation  to  an  electric  dipole,  and  an  electrically  small  loop  is  a good 


approximation  to  a magnetic  dipole.  A third  idealized  probe  of  interest  is  a Huygens'  source  [9] 
which  responds  to  perpendicular  components  of  the  electric  and  magnetic  field.  A Huygens' 
probe  can  be  realized  approximately  by  a small  aperture  antenna,  such  as  an  open-ended 
waveguide  [3],  or  by  a double-gap  loop  [10]  which  simultaneously  measures  perpendicular 
components  of  the  electric  and  magnetic  fields.  This  type  of  idealized  probe  is  closely  related  to 
the  uniqueness  vector  which  is  defined  as  [1 1] 


U = - —n  ^ H, 

tan  ’ 


(16) 


where  » is  the  unit  normal  directed  into  the  region  in  which  the  fields  are  being  expanded.  For  the 
exterior  field  expansion  in  this  section,  n = f.  Specification  of  U over  a closed  surface  uniquely 
determines  the  fields  in  an  adjoining,  source-free  region  (either  interior  or  exterior).  In  theory,  the 
constant  l/r)  can  be  any  nonzero,  real  constant,  but  in  practice  a choice  of  1/t|  (=  \[yJe)  usually 
yields  good  numerical  results  [11].  That  choice  also  corresponds  to  the  response  of  an  idealized 
Huygens'  probe. 

The  spherical  harmonic  expansion  for  U\s  obtained  by  substituting  eqs  (1)  and  (2)  into  eq 
(16)  and  making  use  of  the  tangential  spherical  harmonics; 
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The  orthogonality  relation  in  eq  (1 1)  can  again  be  used  to  solve  for  the  coefficients: 
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It  is  useful  to  examine  the  as5mptotic  form  of  the  square  bracket  factor  in  the  denominator  of  eq 
(18)  to  confirm  that  it  is  well  behaved  numerically.  For  kr  » N , the  radial  functions  can  be 
approximated  [1] 
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Hence  for  large  kr  the  terms  in  the  square  bracket  of  eq  (18)  add  to  give  the  result 
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In  summary,  eqs  (14),  (15),  and  (17)  give  three  methods  of  determining  the  fields  outside 
a sphere  when  the  sources  are  located  inside.  All  three  equations  are  well  conditioned,  and  this  is 
a result  of  uniqueness  properties  that  either  the  tangential  electric  field,  tangential  magnetic  field, 
or  the  uniqueness  vector  U is  sufficient  to  determine  the  external  fields. 

2.2  Exterior  Sources  Only 

This  section  treats  the  special  case  shown  in  figure  2 where  the  sources  are  confined  to  the 
exterior  region,  r > r^.  In  the  source-free  region,  r ^ r^,  the  electric  and  magnetic  fields  can  be 
written  in  terms  of  the  following  expansions  of  standing  waves  that  are  finite  at  the  origin: 

(21) 


(22) 
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For  c = 1,  the  spherical  Bessel  function  which  is  finite  at  the  origin  is  [8] 


(23) 


The  methods  for  determining  the  unknown  coefficients  follow  those  of  the  previous 
section.  Consider  first  the  case  where  the  tangential  electric  field  E ^ specified  over  a surface 
of  radius  r (<  r).  If  a tangential  spherical  harmonic  is  dotted  into  eq  (17),  the  orthogonality 
relationship  in  eq  (1 1)  can  be  used  to  obtain  an  expression  which  is  similar  to  eq  (14): 
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The  only  difference  between  eqs  (14)  and  (24)  is  that  the  outgoing  radial  function  R^\kr)  is 
replaced  by  the  standing  wave  radial  function  R^\kr) . The  derivation  of  eq  (24)  requires  a 
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division  by  and  is  not  valid  for  values  of  s,  n,  and  kr  where  R^\kr)  = 0.  A similar 

procedure  can  be  carried  out  by  starting  with  eq  (22)  for  the  tangential  magnetic  field,  and  the 
result  for  the  coefficients  in  terms  of  the  tangential  magnetic  field  is 


/(-I)"' 


(25) 


This  equation  is  not  valid  for  values  of  5,  n,  and  kr  where  = 0.  These  difficulties  'with  eqs 

(24)  and  (25)  are  to  be  expected  because  either  the  tangential  electric  or  magnetic  field  by  itself  is 
not  sufficient  to  uniquely  determine  the  fields  in  an  interior  region  [11], 

If  both  the  tangential  electric  and  magnetic  fields  are  specified,  then  the  fields  in  an  interior 
region  are  uniquely  determined  [12],  In  such  a case,  eqs  (24)  and  (25)  could  be  used  to  determine 
the  coefficients  because  the  zeros  of  R^^ikr)  and  do  not  occur  simultaneously  for  the 

same  values  of  5,  «,  and  kr.  The  coefficients  can  also  be  determined  from  the  uniqueness  vector 
U as  defined  by  eq  (16).  For  determining  the  fields  in  the  interior  region,  the  unit  normal  in  eq 
(16)  is  given  by  « = - r,  and  the  Huygens'  probe  for  determining  U would  need  to  be  oriented 
appropriately.  The  spherical  harmonic  expansion  for  U is  obtained  by  substituting  eqs  (21)  and 
(22)  into  eq  (16)  and  making  use  of  the  tangential  harmonics: 
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An  equivalent  form  of  eq  (26)  has  been  presented  in  reference  [3],  The  orthogonality  relationship 
in  eq  (1 1)  can  again  be  used  to  solve  for  the  coefficients; 


v/n(-i)" 

,'-'«r[«®(fo-)./*"(fo-)] 


//  £/(r,e,4.)-F  .^(e,<l>)dQ. 


(27) 


Equation  (27)  does  not  suffer  from  the  problem  encountered  in  eqs  (24)  and  (25)  because  the 
denominator  does  not  have  any  zeros  for  real  values  of  kr.  Practical  consequences  are  that  ideal 
dipole  probes  are  not  good  for  determining  interior  fields  and  probes  that  are  similar  to  a 
Huygens'  probe  (respond  to  both  tangential  electric  and  magnetic  fields)  are  preferred. 

2.3  Interior  and  Exterior  Sources 

This  section  treats  the  most  general  case  shown  in  figure  3 where  sources  are  present  in 
both  the  interior  region,  r < and  the  exterior  region,  r > r^.  In  the  source-free  region. 
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r.  ^ r the  electric  and  magnetic  fields  can  be  written  as  expansions  of  standing  waves  and 
outgoing  waves: 


(28) 
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Equations  (28)  and  (29)  can  be  thought  of  as  a superposition  of  eqs  (1)  and  (21)  and  eqs  (2)  and 
(22).  The  coefficients  of  the  standing  waves  characterize  the  exterior  sources,  and  the 

coefficients  of  the  outgoing  waves  characterize  the  interior  sources. 

Since  twice  as  many  unknowns  need  to  be  determined,  more  information  is  needed.  If 
both  the  tangential  electric  and  magnetic  fields  are  specified  over  a sphere  of  radius  r in  the 
source-free  region,  the  orthogonality  relationship  in  eq  (1 1)  can  be  used  to  derive  the  following 
simultaneous  equations: 


(30) 


4n 


The  solution  of  eqs  (30)  and  (31)  can  be  simplified  by  using  the  Wronskian  for  the  radial  functions 


[1]: 


(32) 


From  eqs  (30)  through  (32),  the  solution  for  the  coefficients  can  be  written 
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The  results  in  eqs  (33)  an(J  (34)  are  applicable  to  the  case  where  two  spherical  scans,  one  with  an 
electric  dipole  probe  and  one  with  a magnetic  dipole  probe,  have  been  performed  at  the  same 
radius. 

The  solution  for  the  coefficients  can  also  be  obtained  in  terms  of  two  uniqueness  vectors, 
and  U^,  defined  as 


v,‘  m^dU,  - E^  - -fxS.  (35) 

n h 

is  obtained  by  setting  n = - f and  is  associated  with  inward  propagating  waves.  U^is 
obtained  by  setting  n = f and  is  associated  with  outward  propagating  waves.  From  eqs  (28), 
(29),  and  (35),  their  spherical  wave  expansions  are 
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The  orthogonality  relationship  in  eq  (1 1)  can  be  applied  to  eqs  (36)  and  (37)  to  derive  the 
following  simultaneous  equations  for  the  coefficients: 
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From  eqs  (32),  (38),  and  (39),  the  solution  for  the  coefficients  can  be  written 
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The  results  in  eqs  (40)  and  (41)  are  applicable  to  the  case  where  two  spherical  scans,  one  with  a 
Huygens'  probe  measuring  j/^and  one  with  an  oppositely  directed  Huygens'  probe  measuring  U^, 
have  been  performed  at  the  same  radius.  For  example,  if  an  open-ended  waveguide  is  used  to 
approximate  a Huygens'  probe,  is  measured  by  an  outward  pointing  waveguide,  and  U^is 
measured  by  an  inward  pointing  waveguide. 

For  the  case  where  the  scanning  sphere  is  electrically  large,  approximations  to  eqs  (40) 
and  (41)  for  large  kr  are  useful  for  providing  some  physical  insight.  For  kr»n,  the  required 
radial  functions  can  be  approximated  by  [1,8] 
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From  eqs  (44)  and  (45),  the  coefficients  of  the  integrals  in  eq  (40)  to  order  kr^  are 
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kr 


(46) 


- 0. 
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Thus  the  coefficient  of  the  U.  integral  is  small,  and  eq  (40)  can  be  approximated  by 
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The  physical  interpretation  of  eq  (48)  is  that  can  be  determined  approximately  from  a single 
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scan  with  a Huygens'  probe  that  responds  to  and  not  to  the  outgoing  waves  as  represented  by 
U^.  However,  it  should  be  kept  in  mind  that  eq  (48)  is  valid  for  kr  » n,  and  a complete  spherical 
scan  normally  requires  that  n values  up  to  somewhat  greater  than  kr.  So  eq  (48)  will  not  be  valid 
for  all  relevant  values  of  n. 

Similarly,  the  approximations  in  eqs  (42)  and  (43)  can  be  used  to  approximate  the 
coefficients  of  the  integrals  in  eq  (41); 


(1 


kr 


(49) 
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With  these  approximations,  eq  (41)  can  be  approximated  by 
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(3)  r 


v/n(-ir 


■n 
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4n 

4n 


(51) 


This  result  is  less  advantageous  than  eq  (48)  because  both  the  17^  and  the  integrals  contribute 
to  the  result.  If  the  field  consists  primarily  of  outgoing  waves  (due  to  interior  sources),  then 
and  its  associated  integral  will  be  small. 


3.  LOW-FREQUENCY  EMITTERS 

The  case  of  electrically  small  emitters  is  of  interest  in  electromagnetic  interference  and  in 
detection  of  electronic  timing  devices  for  bomb  detonators.  In  such  cases  the  source  that  is  to  be 
characterized  or  detected  can  be  represented  by  electric  and  magnetic  dipole  moments,  and  the 
characterization  of  the  source  involves  the  determination  of  three  orthogonal  components  of  the 
electric  dipole  moment  and  three  orthogonal  components  of  the  magnetic  dipole  moment. 
Methods  of  determining  these  dipole  moments  using  a three-loop  device  [2]  or  a transverse 
electromagnetic  (TEM)  cell  [13]  have  been  developed,  but  here  the  possibility  of  using  field 
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probes  will  be  analyzed. 

3.1  Electric  and  Magnetic  Dipole  Fields 

Consider  the  geometry  in  figure  3 that  was  analyzed  in  section  2.3  for  general  sources. 

For  electrically  small  sources  and  scan  radii,  the  following  condition  holds;  kr  ^ \ . Since 
r > r.,  kr^is  also  small.  Under  these  conditions,  the  n summations  can  be  truncated  at  « = 1,  and 
eqs  (28)  and  (29)  reduce  to 


2 1 


f-l  m.-\ 


(52) 


2 1 


H(.r,  E 0,  <!>) . 0,  <i>)  ] . (53) 


f-l  TO--1 


The  coefficients  characterize  the  electrically  small  source  of  interest,  and  the  coefficients 
characterize  the  interference  due  to  external  sources. 

The  fields  of  an  electrically  small  source  can  be  considered  to  be  those  of  equivalent 
electric  and  magnetic  dipole  moments,  and  d^,  located  at  the  origin  as  indicated  in  figure  4. 
Each  dipole  moment  can  be  written  in  terms  of  its  rectangular  components: 

d^dx^dy^dimdid  = d x^d  y *d  z.  (54) 

€ «r  9)r  dz  m me  my'  mz  ^ ' 


The  spherical  mode  coefficients  can  be  written  in  terms  of  the  dipole  components  as  follows  [1]: 
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Thus  the  unknowns  of  interest  are  either  six  mode  coefficients  or  six  dipole  components.  The 
radiated  power  P can  be  written  in  terms  of  either  the  spherical  dipole  mode  coefficients  or  the 
dipole  components; 
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(56) 


2 f-i  m—i  2 Stix\ 


3.2  Dipole  Components 

Since  the  number  of  unknown  dipole  components  (or  outgoing  spherical  modes)  is  only 
six,  a dense  spherical  scan  is  not  needed.  By  selecting  the  locations  and  field  components 
properly,  the  equations  for  the  unknown  dipole  components  and  the  unknown  coefficients  of  the 
standing  wave  modes  decouple  to  allow  simple  solutions.  For  electrically  small  sources,  both  the 
radial  and  tangential  fields  are  useful  because  they  are  of  the  same  order  of  magnitude. 

3.2. 1 Solution  in  Terms  of  Radial  Fields 

The  radial  components  of  the  electric  and  magnetic  fields  can  be  obtained  by  taking  the 
dot  product  of  eqs  (52)  and  (53); 


£ (r,0,<|))  . Y, 


A 


m-1 


(57) 


(58) 


Equations  (57)  and  (58)  are  decoupled  because  eq  (57)  includes  only  electric  (TM)  modes  and  eq 
(58)  includes  only  magnetic  (TE)  modes.  A simple  method  for  the  solution  of  eqs  (57)  and  (58) 
makes  use  of  the  values  of  the  spherical  functions  on  the  positive  z axis  (9  = 0 ): 
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Since  four  of  the  six  functions  are  zero,  both  eqs  (57)  and  (58)  can  be  reduced  to  three  two-by- 
two  equations. 

Consider  first  the  radial  electric  field  in  eq  (57).  Only  the  z-directed  electric  dipole 
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produces  a radial  electric  field  on  the  z axis  as  shown  in  figure  5.  If  the  radial  electric  field  is 
measured  at  two  radii,  r,  and  Vj,  then  eq  (57)  yields  the  following  two  equations  in  two 
unknowns: 
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Simultaneous  solution  of  eqs  (60)  and  (61)  yields 
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represents  the  external  sources  and  is  not  of  direct  interest  except  to  indicate  the 
strength  of  the  interfering  fields.  ® is  proportional  to  the  z-directed  dipole  moment,  and  from 

eqs  (55)  and  (63)  is  given  by 
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Similarly,  can  be  determined  from  two  measurements  of  the  radial  electric  field  on  the  x axis, 
and  d^  can  be  determined  from  two  measurements  of  the  radial  electric  field  on  the  axis: 
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The  determination  of  the  magnetic  dipole  components  follows  a similar  method  starting 
with  eq  (58)  for  the  radial  magnetic  field.  From  two  measurements  of  the  radial  magnetic  field  on 
the  z axis,  the  following  two  coefficients  can  be  determined: 
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where  represents  the  external  sources,  and  ig®  is  proportional  the  z-directed  magnetic 
dipole  component  From  eqs  (55)  and  (69),  is  given  by 
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Similarly,  d^  and  d can  each  be  determined  from  two  measurements  of  the  radial  magnetic  field: 
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For  At  « 1 , the  spherical  Bessel  functions  can  be  replaced  by  their  small-argument 
approximations  [8]; 
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For  small  kr^  and  kr2,  eq  (73)  can  be  used  to  reduce  the  dipole  moment  expressions  to 
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To  illustrate  the  immunity  of  this  dipole  formulation  to  interfering  external  sources, 
consider  a plane-wave  field  (from  a distant  external  source)  propagating  through  the  source-free 
region  where  the  radial  electric  and  magnetic  fields  are  being  measured.  (A  single  plane  wave  is 
not  necessarily  the  best  example  of  an  interference  field,  but  it  is  easy  to  analyze.)  For 
propagation  in  the  positive  x direction  and  electric  field  polarization  in  the  direction,  the  electric 
field  Ep  and  the  magnetic  field  Hp  are  given  by: 
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and 


If  these  plane-wave  fields  are  substituted  into  the  dipole  moment  solutions  in  eqs  (74)  through 
(79),  all  of  the  results  are  zero:  d ^ d = d = d = d = d = 0.  The  form  of  eqs  (74) 
through  (79)  is  essentially  that  of  a gradiometer  antenna  which  responds  to  the  gradient  (or 
difference)  of  the  field  rather  than  the  field  itself  This  property  allows  rejection  of  interference 
distant  sources  and  simultaneous  detection  of  nearby  sources. 

3.2.2  Solution  in  Terms  of  Tangential  Fields 

An  advantage  of  using  the  tangential  (transverse)  fields  is  that  measurements  of  the 
electric  and  magnetic  fields  taken  at  a single  radius  r are  sufficient  to  determine  the  12  unknown 
coefficients  in  eqs  (52)  and  (53).  By  using  various  symmetries,  it  is  possible  to  determine  two 
dipole  moments  (out  of  six)  at  a time.  For  example,  consider  the  determination  of  and 
from  two  electric-field  and  two  magnetic-field  measurements  as  shown  in  figure  6.  From  eqs  (3), 
(5),  (52),  and  (53),  the  relevant  field  components  can  be  written: 
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The  solutions  for  the  coefficients  of  the  outgoing  waves  are  obtained  by  taking  sums  and 
differences  of  the  preceding  equations  in  much  the  same  as  is  done  with  the  three-loop  sensor  [2]. 
Using  the  Wronskian  in  eq  (32),  the  results  are 
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Similar  solutions  could  be  obtained  for  the  coefficients  of  the  standing  waves  (due  to  external 
sources).  The  dipole  moments  can  be  written  in  terms  of  the  outgoing  wave  coefficients  by  using 
eq  (55); 
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From  eqs  (85)  through  (87),  the  dipole  moments  can  be  written  in  terms  of  the  measured  fields; 


J\ikr) 


2 2 


\ 


/nifiiV) 


/ \ 
’2  2 


/ J 


(88) 


d - ilTir^ 
ty 


/ ^ 

£.(r,io)-£,(r,i  71) 
2 2 


0)-£,(r,i  n) 
2 2 


\ 


/ J 


(89) 


The  four  remaining  dipole  components  can  be  obtained  from  eqs  (88)  and  (89)  by 
rotation; 
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For  the  simpler  case  where  there  are  no  external  sources,  the  determination  of  the  dipole  moments 
does  not  require  both  electric  and  magnetic  fields.  For  example,  the  expressions  for  the  dipole 
moments  in  terms  of  electric  field  values  are  given  in  the  Appendix. 

For  small  values  of  kr,  the  radial  functions  in  eqs  (88)  through  (93)  can  be  replaced  by 
their  small-argument  approximations  given  in  eq  (73).  Then  the  dipole  expressions  reduce  to 
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To  illustrate  the  immunity  of  this  dipole  formulation  to  external  sources,  consider  again 
the  plane-wave  incident  field  (from  a distant  external  source)  as  given  by  eq  (80).  If  eq  (80)  is 
substituted  into  eqs  (94)  through  (99),  the  following  results  are  obtained: 
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Four  of  the  dipole  components  are  identically  zero  because  of  the  polarization  of  the  incident 
wave.  Both  and  d^^  are  proportional  to  kr  which  is  assumed  to  be  small.  It  can  be  shown  that 
for  arbitrary  incidence  angle  and  polarization,  all  of  the  dipole  components  are  proportional  to  kr. 


4.  CONCLUSIONS 

Spherical  scanning  formulations  have  been  presented  for  the  general  case  where  sources 
are  present  both  inside  and  outside  the  measurement  sphere.  The  spherical  wave  representation  in 
the  source-free  measurement  region  includes  both  outgoing  waves  (due  to  interior  sources)  and 
standing  waves  (due  to  exterior  sources).  The  coefficients  of  both  sets  of  waves  have  been 
derived  in  terms  of  the  transverse  electric  and  magnetic  fields.  Related  results  have  also  been 
derived  in  terms  of  a uniqueness  vector  (a  particular  combination  of  perpendicular  electric  and 
magnetic  fields)  which  can  be  sensed  with  an  ideal  Huygens'  probe. 

The  special  case  of  a low-frequency  emitter  (such  as  a timing  device  for  explosives)  has 
been  studied  in  greater  detail.  The  spherical  wave  expansion  of  this  electrically  small  source  can 
be  truncated  at  « = 1,  and  the  unknowns  are  three  scalar  components  of  the  effective  electric 
dipole  moment  and  three  scalar  components  of  the  effective  magnetic  dipole  moment.  Six 
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additional  coefficients  of  standing  waves  (for  « = 1)  represent  the  interfering  fields  of  external 
sources.  Formulations  for  the  twelve  unknowns  are  derived  in  terms  of  twelve  values  of  either  the 
radial  or  tangential  fields.  The  example  of  an  interfering  plane  wave  (due  to  a distant  source) 
shows  that  the  formulations  are  successful  in  nulling  out  interference  fields. 

A number  of  extensions  to  this  work  are  possible.  The  solutions  in  terms  of  transverse 
electric  and  magnetic  fields  or  uniqueness  vectors  are  representative  of  ideal  electric  and  magnetic 
dipole  probes  or  an  ideal  Huygens'  probe.  A more  general  and  more  practical  measurement 
method  would  treat  multiple  scans  with  a real  probe  antenna  and  its  associated  probe  correction 
[1,3].  The  multiple  scans  could  be  carried  out  at  different  sphere  radii  or  with  probes  pointing  in 
and  out.  Sensitivity  is  a major  issue  for  detection  and  characterization  of  weak,  low-frequency 
emitters.  Calculations  and  measurements  would  be  useful  for  determining  practical  limitations  for 
real  sources  in  a noisy  environment. 
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Appendix:  Dipole  Components  in  Terms  of  the  Transverse  Electric  Field 

For  the  simpler  case  where  no  standing  waves  (due  to  exterior  sources)  exist,  there  are 
only  six  unknown  dipole  components  to  determine.  Following  the  method  of  Section  3.2.2,  the 
six  transverse  electric  components  provide  sufficient  information  for  a solution.  Consider  eqs 
(80)  and  (81)  for  the  special  case  where  the  standing  wave  coefficients,  q!^\  ,,  and  Q^l], 
are  zero.  By  taking  the  sum  and  difference  of  the  simplified  equations,  the  following  solution  for 
the  coefficients  of  the  outgoing  waves  is  obtained: 
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The  corresponding  dipole  components  can  be  obtained  by  substituting  eqs  (Al)  and  (A2)  into  eq 
(55): 
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The  remaining  four  dipole  components  can  be  obtained  from  eqs  (A3)  and  (A4)  by 
rotation: 
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Similar  equations  could  be  obtained  for  the  dipole  components  in  terms  of  magnetic  field  values 
by  starting  with  eqs  (83)  and  (84). 

Equations  (A3)  through  (A8)  are  simpler  than  eqs  (88)  through  (93),  and  they  offer  the 
measurement  advantage  of  requiring  only  electric  field  samples.  They  do  not  reject  the  w = 1 
standing  waves,  but  by  going  through  a plane-wave  example  as  in  Section  3.2.2,  they  can  be 
shown  to  offer  plane-wave  rejection  comparable  to  eqs  (88)  through  (93).  The  tradeoffs  between 
the  two  formulations  merit  further  study. 
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Figure  1 . Geometry  for  a measurement  sphere  of  radius  r with  sources  located  inside  a sphere  of 
radius  r^. 


Figure  2.  Geometry  for  a measurement  sphere  of  radius  r with  sources  located  outside  a sphere 
of  radius  r^. 
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Figure  3.  Geometry  for  a measurement  sphere  of  radius  r with  sources  located  inside  a sphere  of 
radius  and  outside  a sphere  of  radius  r^. 


Figure  4.  Geometry  for  a measurement  sphere  of  radius  r with  equivalent  electric  and  magnetic 
dipoles  located  inside  a sphere  of  radius  r,  and  exterior  sources  located  outside  a sphere 
of  radius  r^. 
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Figure  5.  Geometry  for  a r-directed  electric  dipole  located  at  the  origin  and  exterior  sources 
located  outside  a sphere  of  radius 


Figure  6.  Electric  and  magnetic  field  samples  for  determining  two  dipole  components  in  the 
presence  of  exterior  sources. 
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